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^ ■ We review recent results on QCD at finite temperature. Main emphasis is put on a 
discussion of observables which are of immediate interest to experimental searches for the 
G> ! Quark Gluon Plasma, i.e. the phase transition temperature, the equation of state in 2 
^ ■ and 3-flavour QCD, the heavy quark free energy and thermal effects on hadron masses. 
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1. Introduction 



' Since many years the thermodynamics of quarks and gluons, the transition between 

the low temperature hadronic phase and the quark-gluon plasma (QGP) as well as non- 
perturbative properties of the high temperature phase, have been analyzed through nu- 
merical calculations within the framework of lattice regularized QCD. In particular, the 
^ I studies performed in the pure gauge sector did provide a rather detailed picture of the high 
temperature plasma phase as a medium of strongly interacting partons which are strongly 
screened and influenced by non-perturbative effects even at fairly high temperatures. 
^ . Over recent years thermodjTiamic calculations on the lattice have steadily been im- 

proved. This partly is due to the much improved computer resources, however equally 
^ ; important has been and still is the development of improved discretization schemes, i.e. im- 
^ ' proved actions. This is of particular relevance for thermodynamic calculations which are 
not only sensitive to the long distance physics at the phase transition but also probe 
properties at short distances in calculations of e.g. the energy density or the heavy quark 
potential. The advantages of improved actions for thermodynamic calculations have been 
discussed at recent lattice conferences [|1],0], this does include improvements of the stan- 
dard Wilson and staggered fermion actions as well as new developments of chiral actions 
like the domain wall and overlap fermion formulations. We thus will not go into details 
of these formulations here. We rather will concentrate on a discussion of results, which 
may be of direct interest to experimental studies of the quark-gluon plasma. We will 
discuss the QCD transition temperature and equation of state and will emphasize the 
quark mass and flavour dependence of these quantities which gives some insight into the 
mechanisms controlling the transition to the Quark Gluon Plasma. Furthermore we will 
discuss thermal modifications of the heavy quark free energy and meson masses. 

We also refrain from discussing details of the ongoing investigations of the order of the 
QCD phase transition 0,^, which is first order in the case of three light, degenerate quark 
flavours and most likely is second order in the case of 2-fiavour QCD, although in this 
latter case questions concerning the universality class of the transition and, in particular. 
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the role of the axial anomaly for the transition are not yet settled completely. Our current 
understanding of the QCD phase diagram of 3-flavour QCD at vanishing baryon number 
density is shown in Figure ||. 
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Figure 1. The QCD phase diagram of 3-flavour QCD with degenerate (u,d)-quark masses 
and a strange quark mass rris- 

2. The Critical Temperature and the QCD Phase Diagram 

The early calculations of the QCD transition temperature, which have been performed 
with standard Wilson |^ and staggered fermion actions, led to significant discrepancies 
of the results. These differences strongly diminished in the newer calculations which are 
based on improved Wilson fermions (Clover action) [^^,|^, domain wall fermions ^ as 
well as improved staggered fermions (p4-action) . A compilation of these newer results 
is shown in Figure § for various values of the quark masses. To compare results obtained 
with different actions the results are presented in terms of a physical observable, i.e. the 
ratio of the lightest pseudo-scalar and vector meson masses {mps/my). In Figure ^ we 
show Tc/rriv obtained for 2-flavour QCD while Figure ^ gives a comparison of results 
obtained with improved staggered fermions [^] for 2 and 3-flavour QCD. Also shown there 
is a result for the case of (2+l)-fiavour QCD, i.e. for two light and one heavier quark 
flavour degree of freedom. Unfortunately the quark masses in this latter case are still too 
large to be compared directly with the situation realized in nature. We note however, 
that the results obtained so far suggest that the transition temperature in (2+l)-flavour 
QCD is close to that of 2-fiavour QCD. The 3-fiavour theory, on the other hand, leads to 
consistently smaller values of the critical temperature, Tc{nf = 2)— Tc(raj = 3) ~ 20 MeV. 
Extrapolation of the transition temperatures to the chiral limit gave 
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Figure 2. Transition temperatures in units of my. In Fig. 2a we give a collection of results 
obtained for 2-fiavour QCD with various fermion actions while Fig. 2b shows a comparison 
of results obtained in 2 and 3-flavour QCD with unimproved and improved staggered 
fermion actions. All results are from simulations on lattices with temporal extent Nr = 4. 
The large dot drawn for nips/mv = indicates the result of chiral extrapolations based 
on calculations with improved Wilson [^] as well as improved staggered fermions. The 
vertical line shows the location of the physical limit, mps = m^, my = mp. 



Here rup has been used to set the scale for Tc. Although the agreement between results 
obtained with Wilson and staggered fermions is striking, one should bear in mind that 
all these results have been obtained on lattice with temporal extent A^^ = 4, i.e. at rather 
large lattice spacing, a ^ 0.3 fm. Moreover, there are uncertainties involved in the ansatz 
used to extrapolate to the chiral limit. We thus estimate that the systematic error on the 
value of Tc/ nip still is of similar magnitude as the purely statistical error quoted above. 

We note from Figure |^ that Tc/my drops with increasing ratio mps/mv, i.e. with 
increasing quark mass. Apparently this does not reflect the general picture we have for 
the quark mass (wi^) dependence of Tc. For instance, a simple percolation picture for the 
QCD transition would suggest that Tc(mg) or better Tc{mps) will increase with increasing 
rriq] with increasing rriq also the hadron masses increase and it becomes more difficult to 
excite the low lying hadronic states. It thus becomes more difficult to create a sufficiently 
high particle/energy density in the hadronic phase that can trigger a phase (percolation) 
transition. Such a picture also follows from chiral model calculations UTOf . 

In order to quantify the quark mass dependence of we ideally should find an ob- 
servable to set the scale for T^, which itself is not dependent on rriq. This clearly is not 
the case for my, which actually vanishes in the limit m^ oo. On the other hand we 
know from studies of the hadron spectrum and string tension in quenched QCD (SU(3) 
gauge theory) that these observables agree with experiment and/or QCD phenomenology 
quite well despite the fact that dynamical light quark contributions are suppressed in this 
limit. It thus may be expected that (partially) quenched hadron masses or the string 
tension do provide an almost quark mass independent scale. In fact, this is what tacitly 
has been assumed when one converts the critical temperature of the SU(3) gauge theory 
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Figure 3. The transition temperature in 2 (filled squares) and 3 (circles) flavour QCD 
versus "mpsl\fo using an improved staggered fermion action (p4-action). Also shown 
are results for 2-flavour QCD obtained with the standard staggered fermion action (open 
squares). The dashed band indicates the uncertainty on T^j \fB in the quenched limit. 
The straight line is the flt given in Eq. 2. 

Tel ^Jo ~ 0.63 into physical unitsQ as has been done also in Figure 1^. 

To quantify the quark mass dependence of the transition temperature one may express 
Tf. in units of ^Jo. This ratio is shown in Figure ^ as a function of mpsj \fo. As can be 
seen the transition temperature starts deviating from the quenched values for mps ^ (6 — 
7)-y/a ~ 2.5 GeV. We also note that the dependence of on mpsj \fo is almost linear in 
the entire mass interval. Such a behaviour might, in fact, be expected for light quarks in 
the vicinity of a 2"'^ order chiral transition where the pseudo-critical temperature depends 
on the mass of the Goldstone-particle like 

T,(m.)-T,(0)~m2//^^ (1) 

For 2-flavour QCD the critical indices are expected to belong to the universality class of 
3-d, 0(4) symmetric spin models and one thus would indeed expect 2//55 = 1.1. However, 
this clearly cannot be the origin of the quasi linear behaviour which is observed for rather 
large hadron masses and seems to be independent of nj. Moreover, unlike in chiral models 
[0 the dependence of on mps turns out to be rather weak. The line shown in Figure ^ 
is a flt to the 3-flavour data, which gave 

(T,/v^), = (r,/v^)o + 0.04(1) X with x = mpsl^f^ . (2) 

For the quark masses currently used in lattice calculations a resonance gas model com- 
bined with a suitably chosen percolation criterion would probably be more appropriate 
to describe the thermodynamics close to T^.. 

^We use here and in the following ^Ja ~ 425 MeV which may be deduced from quenched spectrum calcu- 
lations {vripl ^ = 1.81 (4) [|l2|). Similar values have been obtained from partially quenched calculations 
in 3-flavour QCD §. 
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Figure 4. The pressure in QCD with nj = 0, 2 and 3 hght quarks as well as two light 
and a heavier (strange) quark. For nf ^ calculations have been performed on a A^,- = 4 
lattice using improved gauge and staggered fermion actions. In the case of the SU(3) pure 
gauge theory the continuum extrapolated result is shown. Arrows indicate the ideal gas 
pressure psB as given in Eq. 3. 

3. The Equation of State 

When discussing the equation of state of QCD, e.g. the temperature dependence of 
the energy density (e) and pressure (p), we should at least distinguish three regimes; the 
high temperature regime (T>1.5Tc), the critical region (T ~ T^) and the low temperature 
regime (T<0.9rc). The calculation of e as well as p on the lattice is most difficult below 
Tc where these observables are exponentially suppressed, which is true even for realistic 
pseudo-scalar meson masses, mps ^ Tc. We therefore will concentrate on a discussion of 
the two former temperature regimes. 

At high temperature we expect that p/T^ and e/T^ will asymptotically approach the 
free gas limit for a gas of gluons and n/ quark flavours, 

!5£ = %I = fl6 + li„,,)!^ . (3) 

V 2 V30 ^ ' 

From calculations in the quenched limit, the pure SU(3) gauge theory, we know that 
the high-T ideal gas limit is reached only very slowly [|^. In fact, for T ~ (2 — 4)Tc 
thermodynamic quantities deviate by about 15% from the limiting ideal gas value. This 
deviation is too big to be understood in terms of ordinary high temperature perturbation 
theory which converges badly at these low temperatures [|13|. However, it is naturally 



accounted for in quasi-particle models |T^ and resummed perturbative calculations |U 

A similar deviation from ideal gas behaviour has now been found in simulations of 
QCD with 2 and 3 degenerate quark flavours as well as in a simulation with two light and 
one heavier quark mass The results of this calculation which has been performed 



with an improved gauge and an improved staggered fermion action (p4-action) are shown 
in Figure ^ As the p4-action is known to lead to much smaller cut-off distortions in 
the high-T limit than the standard actions, it also becomes meaningful to compare these 



results with continuum models and perturbative calculations at high temperature [15 



We stress, however, that a flnal extrapolation to the continuum limit still has to be done 
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Figure 5. The energy density in QCD. The left (right) figure shows results from a calcula- 
tion with improved staggered (Wilson) fermions on lattices with temporal extent A^.^ = 4 
{Nr = 4, 6). Arrows in the left figure show the ideal gas values ess as given by Eq. 3. 



for QCD with light quarks. From an analysis of the cut-off dependence of the p4-action 
and the experience gained in the pure gauge sector one expects that the results shown in 
Figure ^ are still systematically below the final continuum extrapolated result. 

The pressure shown in Figure ^ for QCD with different number of flavours as well as 
for the pure SU(3) gauge theory clearly reflects the strong change in the number of degrees 
of freedom in the high temperature phase. Moreover, the dependence of Tc on the number 
of partonic degrees of freedom is clearly visible. In view of this it indeed is striking that 
p/PsB is almost flavour independent when plotted in units of T/T^ (Figure ^). 

Unfortunately, Wilson actions with similarly good high temperature behaviour have not 
been constructed so far. The Clover action does not improve the ideal gas behaviour, i.e. it 
has the same inflnite temperature limit as the Wilson action. Consequently one observes 
an overshooting of the ideal gas limit at high temperature which reflects the cut-off effects 



in the unimproved fermion sector [|T6| . These cut-off effects are, however, unimportant in 
the vicinity of the phase transition where correlation lengths become large. It thus makes 
sense to compare results obtained with different actions in this regime. In Figure ^ we 
show recent results for the energy density obtained with improved staggered^ and Wilson 
fermions. We note that these calculations yield consistent estimates for the energy 



density at Tc 

e,^(6±2)T,^ . (4) 
This estimate also is consistent with results obtained for the energy density from calcu- 



lations with a standard staggered fermion action |17 . 



■^This figure for staggered fermions is based on data from Ref ||Tl[. Here a contribution to e/T* which is 
proportional to the bare quark mass and vanishes in the chiral hmit is not taken into account. 
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4. Heavy quark free energy and screening below Tc 

For the discussion of medium effects on heavy quark bound states in the QGP it is 
of great importance to calculate thermal modifications of the heavy quark colour singlet 
potential. Unfortunately this is not directly accessible to lattice calculations, in particular 
not as a gauge invariant observable. At finite temperature one generally calculates the 
heavy quark free energy from correlation functions of Polyakov loops L(n), 



e 



-V{r,T)/T 



(TrL(0)TrLl"(n)) , rT = \n\N^ , (5) 



where Nt denotes the temporal extent of the lattice. At finite temperature the static 
gg-pair can be in a singlet or octet state and the above correlation function thus only 
provides a thermal average over singlet and octet free energies, 

Although Vi = —8 Vs in leading order perturbation theory, we do in general have to deal 
with both contributions to V{r,T) and cannot directly extract the singlet component. 
Nonetheless the latter will dominate at short distances as well as at low temperature. 

The colour averaged heavy quark free energy defined by Eq. ^ thus yields only indirect 
information for the discussion of thermal effects on heavy quark bound states. Nonetheless 
it refiects basic properties like the strong Debye screening of the gg-pair above T^. This 
has been analyzed in detail in the pure SU(3) gauge theory [|19|. Additional screening 



effects do arise in QCD from the presence of light quarks. In addition to quantitative 
changes in the values of screening lengths the most significant changes occur in the long 
distance behaviour of V{r,T). While in the pure gauge sector screening of the heavy 
quark free energy sets in only above T^, i.e. lim V{r,T) remains finite only for T > T^, 
this is the case at all temperatures as soon as rrig < oo. The temperature dependence of 
this screening effect, which arises from the spontaneous creation of gg-pairs in the heat 
bath, is shown in Figure ^. 

As can bee seen the free energy seems to reach a limiting form for T<0.6 Tc. For 
these temperatures it agrees with the confining Cornell type potential up to distances 
r^/o^ ~ 1.5, i.e. r ~ 0.8 fm. With increasing temperature this point is shifted to smaller 
distances. For T ~ 0.95 Tc screening sets in already for r^/a ~ 0.8 or r ~ 0.4 fm. This 
clearly will also have an impact on heavy quark bound states. However, to quantify this we 
need further information on the short distance properties of the heavy quark free energy 
which should allow to disentangle contributions from the singlet and octet components. 

5. Thermal Masses 

The analysis of the temperature dependence of hadron properties, e.g. their masses 
and widths, is of central importance for our understanding of deconfinement and chiral 
symmetry restoration at Tc. Thermal modifications of the heavy quark potential which 
have been discussed above infiuence the spectrum of heavy quark bound states. Their 
experimentally observed suppression thus is expected to be closely linked to the 
deconfining properties of QCD above Tc ||2l|. Changes in the chiral condensate, on the 
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Figure 6. Temperature dependence of the heavy quark free energy in 3-fiavour QCD 



1^. The band of sohd curves shows the Cornell potential V{r) = —a/r + ar with 
a = 0.25 ± 0.05. The finite temperature free energies have been normalized to this 
potential at the shortest distance available, i.e. at rT = 0.25. 



other hand, influence the light hadron spectrum and may result in experimental signatures, 
for instance in the enhanced dilepton production observed in heavy ion experiments . 

In numerical calculations on Euclidean lattices one has access to thermal Green's func- 
tions Gh{t, r) in fixed quantum number channels H, to which in particular at high tem- 
perature many excited states contribute. As the temporal direction of the Euclidean 
lattice is rather short at finite temperature one usually has restricted numerical investi- 
gations to the analysis of the long distance behaviour of spatial correlations functions, 
GH{T,r) ~ exp(— |r|), which defines hadronic screening masses rfiH- This indeed 
gives evidence for the restoration of chiral symmetry above T^, e.g. one finds that scalar 
and pseudo-vector screening lengths become degenerate and also the difference between 
screening lengths in scalar and pseudo-scalar channels strongly diminishes, which gives 
indications for a partial restoration of the Ua{^) symmetry 

In order to get information on the T-dependence of pole masses and their widths one 
has to analyze the structure of temporal correlation functions. The information on hadron 
masses and quasi-particle excitations is then encoded in the spectral function cr/f(u;,p), 

Gnir^p) = [ exp(zpf) GH(r,f) = H duj auiu^.P) '"'^.^'t^,'' ~ ■ (7) 

J Jo smh(u;p/2) 

At finite temperature the temporal correlation function usually is determined only at a 
small number of lattice grid points as the temperature is related to the finite extent of the 
lattice in this direction, A^^ = l/(aT). A way out may be the use of anisotropic lattices 
These calculations indeed show large changes in meson correlators in the vicinity 



of Tc. To what extent the results suggest that pole masses in light meson channels exist 
as well defined states even above Tc, however, is difficult to judge solely on the basis of 
standard analysis techniques also used for correlation functions at zero temperature. At 
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Figure 7. A first attempt to reconstruct the thermal pseudo-scalar spectral function (right) 
from the temporal correlation function (left) calculated within the static approximation 
of QCD at T = O.eTc and ST^ on lattices of size 24^ x 16 and 64^ x 16, respectively p8 . 



high temperature refined statistical techniques like the maximum entropy method pj] may 
be of help. It has recently been shown that this does allow to extract directly the spectral 
functions of hadron correlators [^,^] at T = as well as high temperature |2^. A first 
attempt to use this approach for the determination of the pseudo-scalar spectral function 
below and above is shown in Figure |^. This clearly refiects the drastic, qualitative 
change which occurs in the pseudo-scalar channel when one crosses T^. The Goldstone 
pole which is present in the spectral function below disappears in the plasma phase. 
Of course, this has to be analyzed more quantitatively in the future. 



6. Conclusions 

We have focused in this review on the calculation of basic thermodynamic quantities 
which are of immediate interest to experimental searches for the Quark Gluon Plasma. 

Calculations of the QCD transition temperature which are based on different discretiza- 
tion schemes do provide a rather consistent picture. Extrapolations to the chiral limit are 
in good agreement and suggest a transition temperature of about 175 MeV in 2-fiavour 
QCD. The numerical results obtained so far suggest a rather weak quark mass and fiavour 
dependence of Tc. This also holds true for the critical energy density, which for 2-fiavour 
QCD is estimated to be ec — 700 MeV and thus is of similar magnitude as that found 
already in the pure gauge sector. However, the error on ec still is about 50% and mainly 
is due to the current uncertainty on which is estimated to be about 10%. 

An analysis of the heavy quark free energy below shows that screening effects in- 
crease significantly already below Tc. While the spontaneous creation of gg-pairs leads 
to screening of the heavy quark potential at r ~ 0.8 fm at low temperature this starts 
already at r ~ 0.4 fm for T ~ 0.95rc. 

We do have reached some understanding of thermal effects on hadron properties. In 
particular, modifications of the light meson spectrum due to fiavour and approximate 
Ua{^) symmetry restoration have been established. Indications for medium effects on 
pole masses have been found in an analysis of spectral functions. However, at present 
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these calculations are not detailed enough to be confronted with experimental data. 

Of course, there are many other important issues which have to be addressed in the 
future. Even at vanishing baryon number density we do not yet have a satisfactory 
understanding of the critical behaviour of 2-flavour QCD in the chiral limit and the 
physically realized situation of QCD with two light, nearly massless quarks and a heavier 
strange quark has barely been analyzed. Moreover, the entire phase diagram at non-zero 
baryon number density is largely unexplored. An interesting phase structure is predicted 
in this case which currently is not accessible to lattice calculations. 

Acknowledgements: I would like to thank N. Christ for communication on recent 
thermodynamic results with domain wall fermions. 



REFERENCES 

1. F. Karsch, Nucl. Phys B (Proc. Suppl.) 83-84 (2000) 14. 

2. S. Ejiri, Nucl. Phys B (Proc. Suppl.) 94 (2001) 19. 

3. C. Bernard et al., Phys. Rev. D60 (2000) 054503 and references therein. 

4. A. Ah Khan et al. (CP-PACS), Phys. Rev. D63 (2001) 034502. 

5. K.M. Bitar et al., Phys. Rev. D43 (1991) 2396. 

6. C. Bernard et al., Phys. Rev. D56 (1997) 5584 and references therein. 

7. R.G. Edwards and U.M. Heller, Phys. Lett. B462 (1999) 132. 

8. N. Christ, private communication. 

9. F. Karsch, A. Peikert and E. Laermann, |hep-lat/0012023| . 

10. H. Meyer- Ortmanns and B.-J. Schaefer, Phys. Rev. D53 (1996) 6586; 

J. Berges, D.U. Jungnickel and C. Wetterich, Phys. Rev. D59 (1999) 034010. 

11. F. Karsch, A. Peikert and E. Laermann, Phys. Lett. B478 (2000) 447. 

12. H. Wittig, Int. J. Mod. Phys. A12 (1997) 4477. 

13. B. Kastening and C. Zhai, Phys. Rev. D52 (1995) 7232. 

14. A. Peshier, B. Kampfer, O.P. Pavlenko, and G. SofT, Phys. Rev. D54 (1996) 2399; 
P. Levai and U. Heinz, Phys. Rev. C57 (1998) 1879; 

15. for a review see J.-P. Blaizot and E. lancu, |hep-ph/0101103 . 



16. A. Ali Khan et al. (CP-PACS), |liep-lat/ 0103028 



17. C. Bernard et al., Phys. Rev. D55 (1997) 6861. 

18. L. McLerran and B. Svetitsky, Phys. Rev. D24 (1981) 450. 

19. O. Kaczmarek et al., Phys. Rev. D62 (2000) 034021. 

20. M.C. Abreu et al. (NA50), Phys. Lett. B477 (2000) 28. 

21. H. Satz, Rept. Prog. Phys. 63 (2000) 1511. 

22. G. Agakishiev et al. (CERES), Phys. Rev. Lett. 75 (1995) 1272 and Phys. Lett. B422 
(1998) 405. 

23. Ph. de Forcrand et al. (QCDTARO), Phys. Rev. D63 (2001) 054501. 

24. M. Jarell and J.E. Gubernatis, Phys. Rep. 269 (1996) 133. 

25. Y. Nakaharai, M. Asakawa and T. Hatsuda, Phys. Rev. D60 (1999) 091503 and Eep 



lat/0011040 



26. see also contribution by T. Hatsuda to this conference. 

27. L Wetzorke and F. Karsch, |hep-lat/ 0008008 . 



28. F. Karsch, E. Laermann, P. Petreczky, S. Stickan and L Wetzorke, work in progress. 



